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Abstract We study self-similarity in one-dimensional probabilistic cellular automata 
(PCA) using the renormalization technique. We introduce a general framework for 
algebraic construction of renormalization groups (RG) on cellular automata and apply 
it to exhaustively search the rule space for automata displaying dynamic criticality. 

Previous studies have shown that there exists several exactly renormalizable de- 
terministic automata. We show that the RG fixed points for such self-similar CA 
are unstable in all directions under renormalization. This implies that the large scale 
structure of self-similar deterministic elementary cellular automata is destroyed by 
any finite error probability. 

As a second result we show that the only non-trivial critical PCA are the different 
versions of the well-studied phenomenon of directed percolation. We discuss how 
the second result supports a conjecture regarding the universality class for dynamic 
criticality defined by directed percolation. 

Keywords Renormalization, cellular automata, self-similarity, universality, directed 
percolation. 



1 Introduction 

Elementary PCA, i.e. (l+l)-dimensional, nearest neighbor two-state probabilistic cel- 
lular automata, are often used to model an interesting class of critical dynamic be- 
havior referred to as directed percolation. It has been conjectured (4j that any process 
which 

- has a continuous phase transition from a fluctuating phase into a unique absorbing 
state; 

- is characterized by a positive one-component order parameter near this transition; 
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- has only short range interactions and 

- has no unconventional attributes such as additional symmetries or quenched ran- 
domness 

belong to the directed percolation universality class. In this paper we introduce a gen- 
eral renormalization framework to explore if there are other PC A that show critical 
behavior, which would contradict the stated conjecture. 

In a recent study, Israeli and Goldenfeld Q explored possible coarse-graining of 
elementary deterministic cellular automata (CA). The coarse-graining procedure con- 
sisted of projections of the state space where configurations in local neighborhoods 
of states was mapped onto a coarse-grained neighborhood, or to be more precise 
a coarse-grained light cone. It was shown that many of the elementary cellular au- 
tomata can be mapped onto each other by this procedure. Especially 21 of the 256 
elementary CA are self-similar and can be mapped onto themselves by an appropriate 
projection. It is interesting to ask in which sense the dynamics of these CA can be 
viewed as critical. In this paper we address this question using real-space renormal- 
ization of PCA, which includes the usual CA as deterministic limits. 

2 Renormalization of probabilistic cellular automata 

We first formulate the dynamics of cellular automata in terms of matrices. We then 
present an algebraic framework for renormalization of probabilistic cellular automata 
on square and diamond lattices. 

2.1 Matrix formulation of cellular automata 

Cellular automata (CA) are discrete deterministic dynamical systems with local in- 
teractions, i.e. a collection of cells {<7,(f )} situated on a lattice together with a local 
update rule. Each cell can be in one of a finite number of states {0, 1, ...,S — 1}. The 
update rule is applied synchronously and independently on every neighborhood of 
cells. The concept is easily generalized to probabilistic cellular automata (PCA) by 
considering non-deterministic update rules. 

For our purposes, it is advantageous to use a spin representation and consider a 
matrix form of the update rule. Each local state, or spin, is represented by a vector s of 
length S corresponding to a column of the identity matrix in R s . Collections of spins 
are represented by tensor products of their constituents, {(7i, 02, ...} <-> (s\,S2,---) = 
si (g> *2 ® 

The size of the neighborhood influencing each spin depends on the dimension 
and topology of the lattice. Let the number of such spins be denoted z- We then 
define the matrix P as the S x S z matrix which corresponds to the local update rule 
for a single cell. For deterministic cellular automata, each entry of P is or 1, but 
in the PCA generalization P is a probability matrix (Pjj G [0, 1] and L,-P,'; = 1). An 
update of a spin is achieved by the matrix product where P operates on the influencing 
neighborhood, S((t + 1) = P ■ ILi'j) s j( t ) where the product should be interpreted in 
terms of tensor products and (i,j) denotes neighboring cells. 
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As an example, consider Wolfram's elementary cellular automata 1111 . They have 
two states which we represent with (1,0) if a = 1 and with (0,1) if O = 0. Fur- 
ther, they are defined on a square lattice so every spin s, has three influencing pre- 
ceding neighbors (s' i _ l ,s' i ,s' i+1 ). This neighborhood is a vector of length 2 3 and P 
is a 2 x 2 3 matrix. As P is a probability matrix, its first row has eight parameters 
Pi = (pi,p2,---p%) G [0, l] 8 and its second consists of the corresponding 1 — For 
example, Pi = (0,1,0,1,1,0,1,0) corresponds to rule 90 in Wolfram's nomenclature 
IPTD as can be seen by successive application to the different spin configurations. 

To coarse-grain an automaton updates of larger neighborhoods must be consid- 
ered. These aggregated update rules also have matrix representations, possible to ex- 
press in terms of P. Let the matrix that updates n neighboring spins from their com- 
bined neighborhood be denoted by P„. As the update rule is applied independently 
and in parallel over the lattice, P n can be constructed in terms of a kind of tensor 
product of P:s. However, the neighborhoods will overlap so the usual tensor product 
cannot be used, easily seen e.g. by considering the dimensions of P P • • • P. 

Consider two adjacent spins on a square lattice. They have three influencing 
neighbors each in their nearest neighbor regions. Only four of these are unique as 
the neighborhoods overlap by two spins. Define a 2-overlapping tensor product for 
2 x 2 3 spin matrices A and B as 

(s' u s' 2 | A& 2 B | S ,Si,S 2 ,S3) = (s[ | A | 50,^1,^2) ' (4 I B I *1>'S2,*3) (1) 

in Dirac's vector notation [9 |. With this notation P? can be expressed in terms of the 
basic update matrix as P 02 P. 

More generally, we define a ^-overlapping tensor product for spin matrices A and 
B as 

(ij , ■•• ) Sj| 1 + mi | A (DkB | so,si , ...,s nz ^. m2 -k) = 

(s±,...,S ni | A | So, ...,5„ 2 ) • (<S Bl _|_i, ■■■i^m+m\ I B I s n 2 -ki ■■■■,Sn<i+m.2-k) (2) 

where A (B) describe the transition of n% (mi) spins into n\ (mi) ones. The matrix 
representation of the update rule for n adjacent spins, P„, on a lattice where adjacent 
spins have k overlapping neighbors is now given by a product of n P:s, P„ = P 0^ 
PQ k --P. 



2.2 Coarse-graining of probabilistic cellular automata 

The matrix representation can be used to calculate both coarse-graining and renor- 
malization of probabilistic cellular automata on square and diamond lattices. We 
denote the matrices defining the dynamics by P D and P° respectively. The coarse- 
graining transformation consists of a projection of a block of cells into a single cell 
in space as well as a stroboscopic coarse-graining in time. This combination is needed 
to keep the structure of the light cone. The projection can be written as a S x S N ma- 
trix n. The goal is to calculate the effective dynamics on the coarse-grained level, 
which also has an S x S z matrix representation P. 



Fig. 1 Left: light cones, defined as the influencing neighborhoods of a given cell, on a square and a 
diamond lattice. Right: blocks used for a coarse-graining of a PCA on a square lattice. 



The transformation has two steps. First we define a cellular automaton on blocks 
of cells on the original lattice. This CA will have an alphabet of size S N and each time 
step will correspond to N time-steps of the original CA. We denote the dynamics on 
the block level by Q D on the square lattice and by Q° on the diamond lattice. This first 
step is straightforward and there are no formal requirements on the dynamics. In the 
second step, the alphabet over blocks of states is projected onto the original one-block 
alphabet, which results in a coarse-grained version on the same form as the original 
CA. For the second step to result in a well defined dynamics, i.e. for the state of the 
automaton at time t + 1 to be independent of the state at time t — 1 given the state 
at time t (the Markov property), there are restrictions on both the projection and the 
dynamics ||5]|6]- These restrictions define the coarse-graining and can be formulated 
in terms of the matrices fl and P. Later when we discuss renormalization the coarse- 
graining restrictions are only approximately fulfilled. 

Due to the 45° angle of the light cones on the lattices, « adjacent spins will have 
77 + 1 influencing neighbors on a diamond lattice and n + 2 on a square lattice. For 
a block size of N, the block CA corresponds to N time steps from 3N onto N spins 
on the original lattice for the square case. Figure Q] illustrates the case N = 2. On 
a diamond lattice the influencing neighborhood consists of 2N spins. The transition 
matrix for the block CA, denoted by Q, can be concisely expressed in terms of the 
overlapping tensor product of equation (0. On the diamond lattice we have 



Q° —Pn' Pn+i ' ' ' P 2N-\ (3) 
p n ■ p n+2 ■ ■ ■ P 3N-2- The description of P$ 



and on the square lattice Q D = P^ -^+2 ' ' ' P 3N-2- ^ e description of P^ in terms of 



P D was given in the previous section as an tensor product with overlap 2. On the dia- 
mond lattice the overlap of neighborhoods of adjacent spins is 1, so the corresponding 
F% are given by P^ = P° 0i P° 0i ■■■!*>. 



SxS s~ S (4) 



S N x S N ^ S N 



For a projection fl to form a well defined coarse-grained dynamics, the diagram 
(|4]i must commute in the sense that a time evolution through Q followed by a projec- 
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tion should give the same result as a projection followed by a time evolution through 
P, see |6) for details. Algebraically this means that there must exist an S xS 2 matrix 
P°, representing the coarse-grained dynamics, such that 

n-Q° = p° -(n®n) (5) 

and n ■ Q D = P D ■ (17 <g> 17 ® TI). If P° exists, then equation [5] can be solved using a 
pseudoinvers^H (IT ® 17) + , 

p° = n-Q°-(n<z>n) + , (6) 

and correspondingly P D = TI ■ Q a ■ (17 (g> 17 (g> I7) + on the square lattice. 



2.3 Renormalization of probabilistic cellular automata 

When applying renormalization theory one is usually much less concerned with using 
exact coarse-grainings. As long as a projection which preserves the symmetries of the 
problem is used one usually gets correct quantitative behavior. With this in mind it 
is natural to use equation (0 as an approximate effective dynamics, ignoring that 
equation (O is usually not fulfilled. The renormalization group transformation is 
then given by equations ([3]) and © as 

a N (p>) = n ■ p^ ■ p^ +l ■ ■ ■ p 2 v i • (n ® n)+ (7) 

on the diamond lattice and 

@ N {p°) = n ■ p° ■ p° +2 ■ ■ -p°_ 2 ■ (rj ® n <g) n)+ (8) 

on the square lattice. 

The pseudo-inverse can be interpreted as a back-projection from the coarse-grained 
description to the corresponding states in the original automaton. If diagram [4] com- 
mutes then any such state will result in the same subsequent dynamics. If it does 
not then the way in which the back-projection is performed affects the results of the 
renormalization. The back-projection thus should be done in a way which is consis- 
tent with the dynamics of the original automaton. This is achieved by weighing with 
the stationary distribution over different states. 

Constructing a diagonal matrix Tr± with the stationary distribution over 3A^ — 2 
spins on the diagonal gives us the weighed renormalization transformation on the 
square lattice 

£ N (p) = n ■ p° ■ p° +2 ■ ■ ■ p? N _ 2 ■ (n ® n ® n ■ d° )+ (9) 

On the diamond lattice the stationary distribution is taken over 2N — 1 spins and the 
resulting equation is 

a N {p) = n ■ i% ■ • • ■ p^ n _ r {n®n-D» eq ) + (io) 

The pseudoinverse used is the Moore-Penrose pseudoinverse for rectangular matrices, i.e. for a n x m 
matrix A with m <n: A + = (A T A)^ i A T . 
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Denote the stationary distribution over M adjacent spins by sf q . Since the sta- 
tionary distribution must be invariant under the time evolution, it must satisfy the 
equation 

with / = 1 on the diamond and / = 2 on the square lattice. This is not a closed equation 
as determination of the stationary distribution for M spins require knowledge of the 
distribution for M + 1 spins. The solution lies in approximating the latter in terms of 
the former. This amounts to disregarding statistical correlations over neighborhoods 
larger than M spins. The natural approximatiord is given by 

(si,...,SM)eq-(s2,—,SM+l)eq 

[Sl,S2,—,SM+l)eq= } ; (12) 

[S2, ■■■ 1 S M )eq 

on the diamond lattice and 



/ \ _ (-Slj— ,SAl)eq- (^2,— ,$M+l)eq- {$3,— ,SM+2)eq 

{Sl,S2,---,SM+2)eq — 7 s } \ 

(S2,-.-,S M )eq- (S3,---,SM+l) 



(13) 

eq 



on the square lattice. To calculate the stationary distribution of for example 3N — 2 
spins to order M we solve equation (fTTT) using equation ( fT3l and then sum over spins 
to arrive at s;^~ 2 . This gives increasingly better approximations for larger M, taking 
larger statistical correlations into account. This defines a hierarchy of approximations, 
where the unweighed inverse projections of equations (|7]i and ([8]) can be seen as 
zeroth order. 

We conclude this section by recalling some of the key points of renormalization 
theory. An introduction to the subject is found in for example Q. As macroscopic 
observations of a system correspond to a severely coarse-grained version of the mi- 
croscopic dynamics the fixed points of the transformation are central. For the cellular 
automata these are the ones described by a matrix P* such that 

<% N (P*)=P* (14) 

By repeated application of the renormalization transformation, each reducing 
length scales of the system, it is easy to see that for example the correlation length at 
a fixed point is either zero or infinite. The latter corresponds to a system at criticality 
with self-similarity at all scales. The former corresponds to bulk phases of the system 
where the dynamics are either pure noise, corresponding to a high-temperature limit, 
or frozen, corresponding to a low-temperature limit. These are known as trivial fixed 
points. 

For systems close to a fixed point the behavior is governed by the eigenvalues 
of the Jacobian of the renormalization transformation as can be seen by linearizing 
around the fixed point, 

P = P*+J(P*)-SP + 0{(8P) 2 ) (15) 



It can be shown that this is the maximum entropy distribution over M + 1 spins such that it reduces to 
the distribution over M spins when summed over the first or last spin. 
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where J(P*) is the Jacobian of M^. Components of 8P which lie along eigenvectors 
of / with corresponding eigenvalues |A| < 1 will shrink. Such directions are called 
irrelevant as they eventually disappear after repeated applications of In the same 
way, directions along which the eigenvalues of J have magnitude larger than one are 
called relevant as such deviations grow exponentially and drive P away from P* . The 
relevant degrees of freedom are associated with control parameters that must be tuned 
to achieve critical behavior. 

The universality observed among widely disparate systems close to criticality 
stems from the fact that it is the same eigenvalues (i.e. the same control parameters) 
which drives the systems near the fixed point, regardless of their respective micro- 
scopic dynamics. These eigenvalues determine the critical exponents of systems near 
their critical point. 

2.4 Searching for renormalizable models 

One of the advantages with cellular automata is their discreteness and simplicity. 
It allows us to explore all possible (non-stochastic) renormalization projections for a 
given block size. We can search the space of PCA for fixed points of each such projec- 
tion and in principle exhaustively enumerate all renormalizable cellular automata in a 
given class. This provides a program for determining the possibility of self-similarity 
in any given class of cellular automata. 

We here restrict our exposition to projections of blocks of size two. This is moti- 
vated by the fact that renormalization theory shows that the actual change in length 
scales does not matter as long as the projection preserves the relevant symmetries 
of the problem. In this setting this means that the only possible difference between 
renormalizations of blocks of size two and larger ones lies in possibilities of new sym- 
metry preservations. For example, a block two projection cannot constitute a proper 
majority rule such as is usually used for renormalization of Ising models. This does 
not seem to be a problem as shown by the example of compact directed percolation. 
This universality class shows up even for projections not respecting its o 1 symme- 
try. We have done some investigations of projections of blocks of size three without 
finding any new results. 

For a given probabilistic cellular automaton we want to find fixed points of the 
corresponding renormalization transformation, P* such that P* = S%n(P*~). Here the 
advantage of our approach is apparent: using the approximation of equation (0 and 
its analog on the square lattice we have polynomial equations in the parameters of P. 
This allows us to quickly search for self-similar automata. 

On a diamond lattice with blocks of size two the problem is a system of four 
polynomial equations in four variables of degree five. We were able to use polynomial 
solving software to completely solve this system for each of the 14 possible non- 
trivial projections of two cells into one. We used the hierarchical scheme introduced 
in section l2T2l to refine the approximation for the resulting fixed points which do not 
fulfill the commuting condition of diagram|4] 

On the square lattice the corresponding system has eight equations of degree 
eight in as many variables. Our software cannot solve this. Instead we perform lo- 
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The Domany-Kinzcl Plane The Compact Directed Percolation Plane 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 2 The renormalization flow for the automata in the diamond lattice for two planes in the hyper cube. 
To the left is the Domany-Kinzel (TJ plane, characterized by p2 = pi and p4 = 0. To the right is the plane 
characterized by pi = 1 — p4 = 1, showing the compact directed percolation line. 



cal searches for fixed points seeded with random initial conditions for each of the 
possible projections. 

The algebraic method of equation (O is fast, but may miss self-similar automata 
which are fixed points only with the weighing on the stationary distribution included. 
Using a fine grid (21 points for each dimension, i.e. 21 4 total), we compute the vector 
field induced by equation dTOb on the diamond lattice and find its zeros for each of 
the 14 non-trivial projections. Figure|2]shows the field for the projection 

n { (0 ' 0) -> (16) 
11 \ (o,i), (i,o)(i, i)-n 

for two planes of the relevant hyper-cube [0, l] 4 . 



3 Existence and stability of self-similar automata 

3.1 Overview of results 

Our search shows three kinds of fixed points of the renormalization scheme. They 
are: 

- Trivial fixed points - e.g. rules in which there is no correlation in time, most of 
them corresponding to high temperature limits. 

- Deterministic CA - rules for which all transition probabilities are or 1 . These 
are shown to be unstable and correspond to low temperature limits. 

- Non-deterministic CA - rules with stochastic time evolution, corresponding to 
finite critical temperatures. All such solutions are shown to belong to the directed 
percolation (DP) universality class or the compact DP class. 
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As shown in previous work Q, there is a fairly large number of deterministic au- 
tomata which coarse-grain themselves. Perhaps a bit surprising, there is no such 
abundance of non-deterministic ones. In fact, our search finds no non-trivial non- 
deterministic CA admitting a renormalization with block size 2, apart from models 
of directed percolation or compact directed percolation. 

It is known |4| that the critical manifold for directed percolation ends in a point at 
which the universality class changes, called compact directed percolation. The change 
of universality class is due to an extra OhI symmetry. Our purely algebraic method 
(equations (|7]i and ©) cannot find this point, but when weighing is included (equa- 
tions (0 and (TTQl i) it is identified correctly. We see that extra symmetry present at this 
point that changes its universality class is not a problem for the two block projections, 
which indicate their generality. 



3.2 Trivial fixed points 

In renormalization theory, trivial or weakly coupled fixed points are points where 
there are no correlations (or trivial correlations in a frozen system). 

One example of such a point in our context is an automaton with constant proba- 
bilities over all preceding configurations, P^ = P^jViJ. For such automata the spins 
behave independently and each take a value k with some probability = p 6 [0, 1]. 
It is not surprising that we find such fixed points for all projections. For projections 
which assigns more configurations to one state than to another, the trivial fixed points 
are p = and p = 1, corresponding to rules and 255 for the elementary automata. 
For the majority rules with ties broken in different ways, all p € [0, 1] gives fixed 
points, while for example the projection 



has fixed points for p £ {0, 1/2,1}. 

More generally, any rule where the probabilities only depend on a single spin 
in the preceding neighborhood is a trivial point as there are no spatial interactions 
in such an automaton. No non-deterministic automata fitting this description not be- 
longing to the specific case handled above are found as fixed points of any of the 
renormalizations. However, a number of deterministic automata does; more on those 
below. 



3.3 Deterministic fixed points 

Coarse-graining of deterministic elementary cellular automata is studied extensively 
in Q and also discussed in (6). Among the central results in Q is a diagram showing 
how different such CA coarse-grain each other. Our method reproduce all 21 fixed 
points of this diagram as expected. We can also analyze these fixed points using 
renormalization theory. 




(17) 
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Fig. 3 The self-similar structures generated by deterministic automata are destroyed by any finite noise. 
Time evolution from a single seed for Wolfram's rule 60 without (left) and with perturbations of 0.4% in 
pi and p3 (middle and right). Cyclic boundary conditions are used. 

First of all we see that most self-similar CA are trivial fixed points in the sense 
described above. The exceptions are rules 60, 90, 102 and their 4> 1 symmetries 
as well as rule 150. These rules correspond to critical systems as they build fractal 
structure at all scales. 

However, these fixed points are degenerate in another sense. Analyzing the Jaco- 
bian of the renormalization transformation at these points, we find that all allowed di- 
rections (directions which does not cause the probabilities to leave the interval [0, 1]) 
are associated with eigenvalues |A, | > 1. This means that these CA are critical only 
in the deterministic limit and at large scales the self-similar structure is destroyed by 
any stochastic component of the dynamics. Figure [3] shows an example of this: the 
time evolution of rule 60 builds a Sierpinski triangle but a small perturbation to any 
of the parameters breaks up the large scale structure. 



3.4 Non-deterministic fixed points 

The only non-deterministic fixed points found by our algebraic method are models 
of directed percolation (DP). On the diamond lattice we find only one, known as 
the Domany-Kinzel automaton f lj. This point exists in three versions on the square 
lattice. Figure 3] shows time evolutions of these automata slightly above criticality 
started from single seeds. 
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Fig. 4 Time evolution from a single seed for the three different models of directed percolation found as 
fixed points of the renormalization on a square lattice. The automata shown are slightly above criticality. 
Cyclic boundary conditions are used. 



The fixed point on the diamond lattice of the unweighed renormalization of equa- 
tion © is given byPf « (0.8389,0.6096,0.6096,0). This is for the projection 

f (0,0), (0,1), (1,0)^0 

n -\ (1,1) -+1 (18) 

The projection can be understood as one that preserves the absorbing state of the per- 
colation process. Running the automaton shows that this zeroth order approximation 
underestimates the critical point. We therefore refine our result by using the hierar- 
chy described in section 12.21 Iterating equation (fTTT) combined with equation (fT2l) 
until convergence and weighing the inverse projection with this distribution gives an 
implicit equation for a fixed point. We regard the fixed point condition as an optimal- 
ity criterion and use a simple direct search method, compass search [7], to minimize 
the error. Truncating the hierarchy at M — 6 we get the value for the critical point 

« (0.9447,0.5923,0.5923,0), which turns out to be in better agreement with the 
correct one as is confirmed by simulations^ 

We note that while the percolation of the Domany-Kinzel automaton has been 
studied using a renormalization scheme similar to ours in e.g. IflOl . the larger gen- 
erality in our scheme allows us to see explicitly that the automaton is only a fixed 
point if there exists a completely absorbing state, as shown both by the roots of the 
polynomials and by an analysis of the eigenvectors at the fixed point. 



3 There are several methods to estimate the critical point for numerical directed percolation models. We 
plot log-log diagrams of the occupancy at time t. For a description and discussion, see (3J. 
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On the square lattice there are three directed percolation fixed points, shown in 
figure |U The non-symmetric fixed points correspond exactly to the Domany-Kinzel 
automaton if the square lattice is tilted 45° to give a diamond lattice. In other words, 
the rules are independent of either the left or right preceding neighbor, 

« (0.8389,0.6096,0.6096,0,0.8389,0.6096,0.6096,0) (19) 

for the right skewed automaton in the zeroth approximation and similarly for the left 
skewed one. The corresponding refined approximations are given by the values given 
above for the diamond lattice. The symmetric fixed point is situated at 

PP w (0.8073,0.7679,0.3957,0.2227,0.7679,0.7206,0.2227,0) (20) 

and a refinement with a hierarchy truncated at M = 6 gives the point 

PP w (0.9807,0.9522,0.1623,0.0711,0.9522,0.9267,0.0711,0) (21) 

We note that the zeroth approximation also has a set of spurious fixed points above 
the critical point. These disappear when higher order approximations are used. 

When we include the weighing systematically in the search for the diamond lat- 
tice we find a line of fixed points belonging to the compact directed percolation class. 
These automata are on the form 

P 1 ° = (l,p,l-p,0) (22) 

for any p G [0, 1]. These possess an extra 0f!l symmetry compared to the directed 
percolation processes, placing them in another universality class. For p = 0.5 the 
automaton describes a annihilating random walk with the boundaries interpreted as 
particles. Tuning the parameter biases the walk in the corresponding direction until, 
at p = 1 or 0, the deterministic limit is achieved and the particles move in straight 
lines. 

Studying the eigenvalues of the transformation around the compact directed per- 
colation points shows that they have no irrelevant directions. They have an exactly 
marginally relevant direction along the line defined by equation (l22l as any move- 
ment along this line will result in a new fixed point being visited. 

The aim of the method presented in this paper is not to calculate the critical expo- 
nents for any particular universality class. For completeness we nevertheless exam- 
ine the possibilities of such a calculation for the directed percolation fixed point. We 
know that the spatial correlation length scales algebraically close to the critical point 
as 4 ~ f~ v± . Through renormalization theory we know that v± = \nb/\nX where X 
is the largest eigenvalue of the Jacobian of the renormalization transformation and b 
is the length scaling (equal to the block size N in our case) (2). 

For the zeroth order approximation we can do the differentiation analytically. For 
higher order approximations we calculate the Jacobian numerically. Table Q] shows 
the obtained values for different orders of the approximation when applied to the 
Domany-Kinzel automaton at the point reported above. They should be compared 
with the value v± = 1.096854(4) reported in |4). Note that already using neighbor- 
hoods of size four comes within 10% of the correct value. However, the convergence 
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Table 1 The largest eigenvalue X\ and corresponding critical exponent for the renormalization trans- 
formation with different orders of approximation M. 



M 




v± 





1.297 


2.667 


4 


1.994 


1.004 


6 


1.986 


1.011 


7 


1.984 


1.011 


S 


1.984 


1.012 


9 


1.983 


1.012 



for larger neighborhoods is extremely slow. This is most likely due to the geomet- 
ric scaling of the correlations close to criticality. This makes the method unsuitable 
for accurate determination of the exponents. Using the same procedure for the fixed 
points on the square lattice gives qualitatively similar results but with even slower 
convergence. This is not surprising since the model in this case is embedded in a 
large space of confounding variables. 



4 Summary and discussion 

In this paper we use renormalization theory to show that the deterministic self-similar 
elementary cellular automata correspond to low temperature limits where any finite 
probability of errors destroys correlations at large enough scale. This is analogous 
to the situation for the one-dimensional Ising model which lacks large scale structure 
for all T > 0. The situation is quite different in for example the two dimensional Ising 
model where the critical point is defined by only two relevant parameters, tempera- 
ture and the external field. Large scale structures are insensitive to other perturbations 
of the microscopic dynamics, corresponding to irrelevant directions of the renormal- 
ization transformation. This is why the Ising model defines a universality class and 
can reproduce the scaling behavior of a large variety of experimental systems near 
criticality. We conclude that the self-similar deterministic CA do not define any in- 
teresting universality classes and that any perturbation of their dynamics destroys 
large scale fractal structure. Our analysis shows that there are limits on large scale 
(space-time) structure generated by elementary cellular automata with finite error 
probability. Lindgren |[8] showed that the space correlations of automata is destroyed 
by uniform noise. Our method generalizes this to space-time correlations and noise 
in any part of the rule of the automaton. 

By systematically exploring both the space of probabilistic cellular automata and 
renormalization transformations we show that the only non-trivial universality class 
among two-state, one-dimensional nearest-neighbor automata is directed percolation. 
This results can be seen as further evidence for the directed percolation conjecture 
discussed in the introduction, i.e. that any critical system describable by two-state 
nearest neighbor PCA belongs to the same universality class as DP. Our results pro- 
vide no definite proof of the conjecture as only deterministic projections with a rela- 
tively limited range have been considered. 
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